Active suppression of dephasing in Josephson-junction qubits 
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Simple majority code correcting k dephasing errors by encoding a qubit of information into 2fc + 1 
physical qubits is studied quantitatively. We derive an equation for quasicontinuous evolution of 
the density matrix of encoded quantum information under the error correction procedure in the 
presence of dephasing noise that in general can be correlated at different qubits. Specific design of 
the Josephson-junction circuit implementing this scheme is suggested. 



Josephson qubits are among the most promising 
devices to implement solid state quantum compu- 
tation JiLH]. Quantum manipulations of individ- 
ual |§,5^]|J^,§,§ and coupled qubits has been 
demonstrated experimentally. At present, probably the 
main obstacle to the development of a larger-scale solid 
state quantum logic circuits is presented by decoherence. 
It is therefore important to develop strategies to mini- 
mize the effects of decoherence on the dynamics of the 
qubit systems. 

Known approaches to reduction of decoherence include 
both error-correction and error-avoiding schemes that ei- 
ther employ symmetries of the qubit-environment inter- 
action to create areas of the Hilbert space not affected by 
decoherence 11; 12] or use rapid random dynamic pertur- 



bations of the system to average out the effects of external 
noise [^,0 . The error-avoiding approaches appear to be 
less suitable for the solid-state qubits. Indeed, noise in 
solid-state systems typically does not have any particular 
symmetry and its correlation time is short, so that the 
application of the control pulses within this time-scale, as 
required by the dynamic averaging schemes, is problem- 
atic. This leaves error-correction as the main strategy for 
suppression of decoherence in solid-state qubits. In this 
work, we suggest an implementation of one of the basic 
error-correction algorithms for the suppression of dephas- 
ing errors (which can be expected to be the dominant 
type of errors in soHd-state circuits - see, e.g., Ref. (H), 
and develop its quantitative description. Our scheme em- 
ploys the Josephson-junction qubits that combine charge 
and flux dynamics [ p"5|jl^ J^, and requires only a small 
number of qubit transformations to operate. 

From the perspective of the general theory of error- 
correction, an interesting feature of the scheme consid- 
ered in this work is the possibility of developing its de- 
tailed quantitative description within the realistic model 
of the qubit-environment interaction and analyzing, for 
instance, the effect of the correlations in the noise act- 
ing on different qubits. While discussions of the error- 
correction rely typically in an essential way on indepen- 
dent noise models, environments of the solid-state qubits 



can be to a large degree correlated because of the finite 
distance between the qubits in a circuit. A clear illustra- 
tion of this is provided by the background charge fluctu- 
ations that are the main source of dephasing in charge 
qubits P|6l |lC| , [l7|| . Long-range nature of the Coulomb in- 
teraction creates noise correlations by coupling the qubits 
to the same charge fluctuators. 

We consider specifically the problem of "quantum 
memory" , when the task is to preserve the stationary 
state of the qubit in the presence of dephasing noise. 
The qubit Hamiltonian contains then only the coupling 
to the environment. Under the assumption that the en- 
vironment has many degrees of freedom each of which is 
only weakly coupled to qubits, it can be modeled as an 
ensemble of harmonic oscillators |jl^,|l^,^,^ (see how- 
ever [|l^ ) , so that the Hamiltonian of the qubit register 
is: 
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where = fc[-^mj"(w)am,Lj + h.c.]. Here we assumed 
several independent ensembles of environmental oscilla- 
tors (numbered by to), as needed to model different pro- 
files of spatial correlations of random forces ^j. The 
index j = 1,2... in (|^) numbers the qubits, and coeffi- 
cients \,n^j(uj) are coupling constants of the qubit j to 
the oscillators of reservoir to in the mode u and cre- 
ation/annihilation operators am,un <An uj- Time evolution 
of the "qubits-|-environment" system is described conve- 
niently in the interaction representation with respect to 
the interaction Hamiltonian of Eq. (|l|). The evolution 
operator U{t) can then be calculated explicitly by sep- 



arating the two non-commuting parts. 
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of the qubit-oscillator coupling, and using the fact that 
their commutator is a c-number: 



L/.(t)=exp{z^ "^-f "^ |^A„,,.0-)P}. 



(2) 
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The first term in U (t) represents fluctuating phases ipj (t) 
of the qubit basis states induced by the environmental 
forces ^j (t): <Pj{t) = ^j{t')dt'. The second term, Ur{t), 
results from the renormalization of the qubit parameters 
by the qubit-environment interaction. To see this more 
explicitly, we note that the sum over frequencies lo in 
this exponent has a natural cut-off at large frequencies 
w ~ T""'^, where Tc is the time scale at which environment 
forces acting on different qubits are correlated. For weak 
decoherence we are interested in the time scales much 
larger that Tc- In this regime, the phase represented by 
Ur{t) is dominated by the term that grows linearly with t, 
and can be viewed as arising from the renormalization of 
the qubit energy. Equation (||) shows that such a renor- 
malization includes then the shift of the total energy of 
the register and creation of the qubit-qubit interaction. 
The total energy shift is irrelevant as long as we consider 
an individual register. Neglecting it, we see that Ur{t) 
results from the Hamiltonian evolution with the Hamil- 
tonian 

Hr^-Y.V,,,ai^^<jP, (3) 

3 -J' 

and Vjj' = 2ReX;™,<^(Amj(w)A;,_j-,(w)/tj), if the sum 
over frequencies lu in this expression is converging at low 
frequencies. The qubit-qubit interaction strength Vjj' is 
non-vanishing only if the same reservoir m couples to 
more than one qubit, so that the reservoir forces at 
different qubits are correlated. 

The time evolution with the Hamiltonian , and more 
generally, the evolution operator [/,. in Eq. (|^) represent 
deterministic part of the qubit evolution induced by the 
qubit-reservoir interaction. As a result, it can in prin- 
ciple be compensated for by adjusting the regular (non- 
dissipative) part of the Hamiltonian of the qubit regis- 
ter. This procedure, however, is impractical even in the 
case of constant Hr, since the interaction constants Vjj' 
are apriori unknown and incommesurate quantities. This 
complexity means that a more appropriate approach is to 
treat the time evolution represented by Ur as dephasing 
despite its deterministic character. 

The time evolution of the density matrix p{t) of the 
qubit register is obtained from Eq. (||) through the re- 
lation p{t) = Trenv{U'i{t)a{0)U{t)}, where a is the to- 
tal density matrix of the "qubits-|-environment" system. 
The environment will dephase the qubits if they are pre- 
pared initially in the state p(0) that is uncorrelated with 
the state of the environment, cr(0) = Penvp{0). Assuming 
that the environment is in thermal equilibrium at tem- 
perature Q, and no error correction procedure is applied, 
we get using the standard property of the Gaussian noise: 

p(i)=exp{-l^(^,(i)^,.(t))(a«-a«)(ar)-ap'))} 
id' 

■U^it)p{0)Urit). (4) 



Here we introduced the convention that the bar over cr^ 
operators means that they act on p from the right. Qual- 
itatively, Eq. (Q) shows that the matrix elements of p that 
are further away from the diagonal in the basis decay 
faster. The diagonal elements (on which az~ o-^ — 0) re- 
main constant. In the case of one physical qubit, Eq. (Q) 
gives p{t) = e-<'^''(*)>(i-'^-*-V(0), i-e., the off-diagonal 
elements of p are suppressed with time as e"^^"^ = 
g-f(t) jf the environment density of states is Ohmic, 
i-e-, Y.,n,w |Am(w)p... = g duJUJe-'^''°... , direct evalua- 
tion for e -CI/to gives: P{t) = 2gln[sinh(7rte)/(7rroe)]. 
At large t, when the random force ^ appears (5-correlated, 
P{t) reduces to P{t) = Tt, where T = iirgQ is the de- 
phasing rate. 

One can reduce the effective dephasing rate by the en- 
coding that corrects the phase errors . Generalized 
to k errors, this encoding is: 

a|0)+/3|l)->a|+J^^)+/3|^_^^). (5) 

2fc+l 2/c+l 

In Equation (^) , a bit of quantum information is encoded 
in the state of the 2fc -f 1 physical qubits, and the |±) 
states of each of these qubits are obtained through the 
Hadamard transform H (the 7r/2-rotation around y axis) 
from the |0, 1) states. All of the az operators in the 
dephasing-induced time evolution (H) are changed by H: 
HazH = ax, so that for the states on the right-hand- 
side of Eq. (H) the dephasing looks like transitions be- 
tween the |±) states of each qubit, and can be directly 
detected by measurements in this basis and corrected by 
applying simple pulses returning the qubit into the ini- 
tial state. The error-detecting measurements, however, 
should not destroy the quantum information encoded in 
the state (^, i.e., they should not distinguish the a and (3 
parts of this state. This condition is not satisfied by mea- 
surements on individual qubits but can be satisfied by 
the measurements on pairs of the nearest-neighbor qubits 
comparing their states. Despite the apparent complexity 
of this scheme, it has quite natural implementation in the 
Josephson-junction qubits - see Fig. 1. 

To describe this process quantitatively we assume that 
its measurement/correction part can be done on the 
time scale that is much shorter than the one set by the 
characteristic dephasing rate F. Different terms in the 
environment-induced evolution of the encoded state, Eq. 

, during the time interval T between the successive 
application of the "measurement -|-correction" operations 
can be conveniently classified by the number of qubits 
fiipped during this time interval. In the relevant regime 
of sufficiently short T: P{T) <C 1, the probability am- 
plitude of these terms decreases rapidly when this num- 
ber increases. If we keep only the terms that fiip up 
to k qubits, we see directly from Eq. (||) that the time 
evolution at this level of accuracy (denoted by Uk{T)) 
preserves the superposition of the a and /3 parts of the 
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encoded state: 

UkiT)[a\(B) + /3|e)] - Y.^a\4'g) + pR\^,)]u, . (6) 
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Here index q runs over 2^*^ different register states ob- 
tained from the state |0) = | + ...+) by flipping up to k 
qubits, Uq are the probability amplitudes of these states, 
R\ipq) denotes the state \iljq) with all 2fc + 1 qubits in- 
verted, and |0) = I — ...—). 

The measurements that compare the qubit states in all 
pairs of the nearest-neighbor qubits do not distinguish 
states \4>q) and R\ipq), and therefore also preserve the 
superposition of the a and (3 terms in Eq. (||) . The 2^*^ 
different outcomes ("equal" or "different") of the 2k such 
measurements distinguish all terms with different q in Eq. 

and enable one to decide what qubits were flipped 
during the time interval T. Application of the correcting 
pulses should then bring the state of the qubit register 
back to its initial form (^) so that the encoded quantum 
state does not change in this approximation. The resid- 
ual evolution of the encoded state is associated with the 
possibility that environment flips more that k different 
qubits; for P{T) - precisely fc -I- 1 qubits. Following 
the same steps as above, we see that when fc-l- 1 qubits are 
flipped, the measurement /correction cycle interchanges 
the a and (3 weights in the encoded state (||). Since the 
probability p of this mistake is small, p ^ 1, the encoded 
state changes substantially only on the time scale larger 
than the period T of one error-correction cycle, and its 
evolution on this scale can be conveniently described by 
the continuous equation for the density matrix p^'^^ in the 
basis of I®) and |0) states. The interchange of the a and 
P terms in (H) leads to the following evolution of p' 



(c). 
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P++ = y(/'-- 
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(7) 



Here 7^ = 2p/T is the effective dephasing rate of the 
encoded quantum information, and the superscript (c) 
indicates that p^'^^ is the reduced density matrix in the 
presence of error correction. Thus, our error-correcting 
procedure replaces the dephasing in the individual phys- 
ical qubit with the dephasing of encoded quantum in- 
formation at a smaller rate. Indeed, if one writes Eqs. 

in the rotated basis |©) ± |0), they exphcitly acquire 
the form characteristic for pure dephasing: constant di- 
agonal elements of the density matrix and decay of the 
off-diagonal elements with the rate 7fc. The dephasing 
rate 7^ can be calculated from the evolution operator 
(^. Now we discuss several important limits. 

For fc 1, when the relevant errors flip 2 out of 3 
qubits, we get: 

j>3' 

The first two terms in this expression represent contri- 
bution to dephasing from noise correlations at different 



qubits, while the last term exists also for uncorrelated 
noise. If the noise is J-correlated in time, 71 reduces to 
71 = Tj:^^^,i2V^j, + T^, + TjT,,/2), where F, is the 
dephasing rate in the jth qubit, and Tjji is introduced 
through 2{(fj{t)ipj>{t)) = Tjjd. 

For k = 2, the dephasing rate of the encoded state is: 



72 



(8) 



j>3'>j" 



Since the effective coupling induced by the environment 
- see Eq. (^), flips the qubits only in pairs, it does not 
contribute to 72. If the dephasing noise is (5-correlated 
in time, its space correlations are non-vanishing only for 
the nearest- neighbor qubits, and the corresponding de- 
phasing rates are the same for all qubits, Eq. ^ gives: 
72 = 5Fr2(f2 + rV2), where f = F^+i. 

If the dephasing forces at different qubits are uncorre- 
lated, the encoded dephasing rate can be easily calculated 
for arbitrary k: 
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J2 PjAT)PjAT)-Pj,UT), (9) 



3l>32>--->Jk+l 



and one sees that 7^ decreases exponentially with the 
"degree of encoding" k. When the probabilities P{T) of 
dephasing errors in individual qubits can be expressed 
through the dephasing rate F, Eq. (P) reduces to 7^ = 
F(Fr)'=(2fc + l)!/(2'=fc!(fc + 1)!), ff F is the same for aU 
qubits. 

Exponential suppression of 7^, with k is limited in the 
scheme considered above by possible imperfections of the 
measurement/correction operations. The most impor- 
tant is direct dephasing of the encoded state by measure- 
ments, which, in contrast to correction steps, need to be 
performed each period T. For example, one of the spe- 
cific non-idealities of measurement detectors that leads 
to direct dephasing of the encoded state is residual lin- 
ear response coefficient of the quadratic detectors needed 
to perform pair-wise comparison of the qubit states - see 
Eq. (|ll| ) below. Linear terms couple the detector directly 
to the |±) states of individual qubits and introduce finite 
phase shifts between them. Since the number of required 
measurements is proportional to fc, the rate of introduced 
dephasing should also be proportional to fc, and can be 
denoted as 7fc. The effect of this dephasing on the evo- 
lution of the encoded quantum information is described 
then by adding the usual dephasing term to the equation 
for the off-diagonal element of the density matrix p^'^^ (|^): 



.(c) Ik, (c) 

p+- ^ y(/'-+ 



P+-) -lkp\'_ 



(10) 



Qualitatively, the two types of dephasing processes in Eq. 
( |lO| ) have similar effect of suppressing the fidelity of the 
encoded state, but depend differently on k. The optimum 
degree of encoding is estimated crudely by minimizing 
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the total dephasing rate: kopt ^ lii(7/r)/lii(Tr). One 
obvious result of this optimization is that for the consid- 
ered scheme of the dephasing suppression to make sense, 
the dephasing introduced by imperfections of the correct- 
ing procedure should be much weaker than the original 
qubit dephasing F. 
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FIG. 1. Schematics of the Josephson-j unction circuits im- 
plementing suppression of dephasing. Crosses denote tunnel 
junctions, the electrodes between them act as charge qubits. 
Monitored currents in the nearest-neighbor loops enclosing 
qubits allow to detect dephasing errors. 

This condition can be satisfied in Josephson-junction 
qubits, where the dynamics of magnetic flux character- 
ized by longer coherence times (at least tens of nanosec- 
onds - see, e.g., [^) can be used to suppress dephasing 
in charge-based qubits. The charge qubits have quite 
short decoherence times, ~ 1 ns limited by the 

background charge fluctuations, but offer some advan- 
tages, e.g., demonstrated simplicity of the qubit-qubit 
coupling jl^ . Therefore it would be of interest to use the 
approach discussed in this work to suppress dephasing 
of charge degrees of freedom with the help of controlled 
flux dynamics. A sketch of the possible set-up achiev- 
ing this is shown in Fig. 1. Its main elements are the 
charge qubits, formed by two small tunnel junctions in 
series, enclosed in small superconducting loops threaded 
by magnetic flux $ equal to half of the magnetic flux 
quantum $0. It can be shown [ p6[ that the current in 
each loop represents the component of the qubit dy- 
namics, and its monitoring measures therefore the qubit 
in the ax basis as needed for detection of the dephasing 
errors. Comparison of the states of the nearest-neighbor 
qubits can be achieved by measuring not directly the cur- 
rents in the loops, but the square of the difference (or of 
the sum) between the currents. Such a quadratic detec- 
tion measures the product operator (Ti cri : 



(a(^)±a 



(J+l)^2 



2(l±aW4^-+i)), 



(11) 



and provides information on whether the states of the 
two qubits are the same or not without measuring them. 
Quadratic detection can be realized by the usual mag- 



netometers but operated at a point where the linear re- 
sponse coefficient vanishes. These measurements, subse- 
quent classical calculations, and application of correction 
pulses, can be done with sufficient frequency by existing 
"SFQ" superconductor electronics compatible with 
the qubits. 

In summary, we suggested a simple scheme of perform- 
ing basic error-correction in Josephson-junction qubits. 
The scheme suppresses dephasing errors and can be an- 
alyzed quantitatively within the realistic model of envi- 
ronment, including the possibility of noise correlations at 
different qubits. If the errors introduced by the correc- 
tion procedure are negligible, the residual dephasing rate 
for the encoded quantum information decreases exponen- 
tially with the degree of encoding. 
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